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Abstract. Two-dimensional radiative transfer is employed to obtain the broad-band infrared spectrum of active 
galaxies. In the models we vary the geometry and size of the obscuring medium, the surface density, the opacity 
and the grain size distribution. Resulting spectral energy distributions are constructed for different orientations of 
the toroid. Colour-colour comparisons with observational data are consistent with previous observations that the 
emission longward of 60 (j,m is produced by star-formation and unrelated to the presence of an obscuring torus. 
We also find that the toroid cannot be flat, but is rather conical or flaring. The density is most likely constant 
with radius, and the size is relatively large with an inner radius around 10 pc. A direct comparison with radio 
galaxy Cygnus A yields a best fit for a conical disk with constant surface density, and a size from 10 to 30 pc, 
assuming the far-infrared emission is due to star-formation in the host galaxy. 

Key words. galaxies:active, galaxies :Seyfert, galaxies:starburst, galaxies individual: Cygnus A, infrared:galaxies, 
infrared:ISM, quasars:general, ISM:dust 



1. Introduction 

Unified models for active galaxies postulate the existence 
of a circumnuclear dust torus. This torus obscures the 
central regions from view at certain solid angles, at wave- 
lengths where the dust is optically thick. Various authors 
have studied the behaviour of a dusty torus in the pres- 
ence of a strong radiation field. Models have been con- 
[ structed for the absorption of the optical to FUV emis- 
sion by the dust, and the re-emission of this energy in 
the infrared JPier fc Krolikll992HGranato fc Danes<Jl994l 
lEfstathiou fc Rowan- Robinsonlll995|) . These models have 
been used to interpret IRAS observations of active galax- 
ies. Their main result is that the spectrum from a dusty 
disk is much broader than a single grey body. However, 
observed SEDs are often still broader than the models, 
and show different behaviour of the 10 /xm silicate feature. 

As discussed in previous work l|van Bemmel fc Barthell 

l200l[) . it is very likely that there are two dust compo- 
nents in active galaxies. One component is a large scale 
dust distribution, heated by stars. The other is the ob- 
scuring torus. In narrow-line radio galaxies and quasars it 
is difficult to disentangle the two components, since their 
broad-band spectra are not well known, due to their dis- 
tance and the sensitivity of previous infrared instruments. 
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Previous model studies of obscuring tori were only con- 
cerned with another subset of the active galaxies, namely 
Seyfert galaxies, which happen to be nearby and therefore 
bright and well observed infrared sources. Moreover, these 
studies focused on reproducing the general width of the 
SED. Furthermore, the emphasis was on fitting a few se- 
lected objects, rather than a complete sample. In previous 
studies, the model geometry has been fixed, and the de- 
pendence of the resulting SED on the geometry and other 
disk parameters is not well understood. Finally, the mod- 
els assume standard Galactic dust, which might not be 
appr opriate for luminous active galaxies. Several s tudies 
(e.g. lLaor fc Draffielll993t iMaiolino etH\ l2001albl ) have 
shown that the lack of a 10 /im silicate emission feature 
and shallow 2200 A absorption can point to dust that is 
dominated by relatively large grains, while Galactic dust 
is dominated by small grains. 

More recent models treat the torus as a clumpy 
medium, which produces a broader infr ared spectrum, 
more consistent with the observations l|Nenkova et alJ 
2002j). However, these authors do not vary the geometry 
of the torus, nor the distribution of the mass. Although it 
is likely that the torus does not have a smooth density dis- 
tribution, it is not yet understood how other parameters 
influence the width of the resulting broad-band spectrum. 
Moreover, they assume that all infrared emission observed 
in active galaxies arises in the torus. 
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Parameter 


Symbol 


Values 


Reference 
model 


Equatorial optical depth at 20 /im 


T~20 


5, 10, 20, 50 


10 


Scale height exponent 


a 


0, 1, 9/7 


1 


Surface density exponent 


P 


0, -1, -2 





Inner radius of the disk (pc) 


ro 


3, 30 


3, 30 


Outer radius of the disk (pc) 


Tout 


30, 300 


30, 300 


Aspect angle (degrees) 


e 


5, 30, 50, 60, 70, 90 


all 


Grain size exponent 


q 


-2, -3.5 


-2 


Scale height at inner radius (pc) 




= r 


3 


Spectral index of input spectrum 


j 


-2 


-2 


100-400 eV luminosity 


Lx 


1O 12 L 


1O 12 L 



Table 1. The different values for the model parameters, and their value in the reference model. The inclination is 
defined with respect to the symmetry axis of the torus. 



Here we will re-examine the infrared SEDs of active 
galaxies from a theoretical point of view. We employ a 
new and very versatile radi ative transfer code RADICAL 
l|Dullemond fc Turollalfe oOO) , and examine a wide param- 
eter space. Specifically, we study the dependence of the 
emission on observing angle, optical depth, disk geometry 
and size, dust mass distribution and dust grain size distri- 
bution. The initial goal will be to constrain the parameters 
of the torus, by comparing the colours to observations, and 
by fitting the infrared broad-band spectrum of Cygnus A. 
As a first order approximation we will assume the torus 
has a smooth distribution of dust. Eventually, these new 
models can be used to disentangle the two dust compo- 
nents in active galaxies, which will allow us to study star- 
formation in active galaxies, and the parameters of the 
obscuring torus in more detail than before. 



2. Model setup 

In our disk model a geometry for the dust, an input spec- 
trum and dust opacities need to be specified. The numer- 
ical method is described below in Sec. 01 We have con- 
structed a grid of models, built by varying different pa- 
rameters. This results in a set of over 150 model SEDs, in 
which all variables are cross-correlated. Below we describe 
the different parameters of the models (see also Table EJ. 

2.1. Disk geometry and density distribution 

The density of the disk in angular direction is a Gaussian, 
with a width of h(r). The mass density for any point in 
the disk is then given by: 



exp 



Mr) 



with z — r*t&n(9) and E(r) the surface density. Due to the 
setup of the model, the surface where T20 equals unity does 
not coincide with the scale height of the disk. For higher 
opacities this surface lies at higher z. Therefore the actual 
opening angle of the torus is smaller for higher opacities, 
although the width of the Gaussian is kept constant. 



The geometry of the disk is thus described by four 
parameters: the scale height h(r), the surface density S(r), 
the inner radius r and the outer radius r out . Each of these 
parameters is varied. For the scale height we assume: 



h(r) =h [- 
and the surface density is defined as: 



E(r) = So 



with ro the inner radius of the disk. 

For hg we assume it equals the inner radius, i.e. h(ro) 
h 1 r 1 . so the disk is thicker if it is further away from the 
emission source. This implies that the opening angles are 
comparable for similar values of t 2 o ■ The resulting opening 
angles are consistent with the size of the broad-line region 
that has to be shielded in a more edge-on orientation, and 
with the observed opening angles of ionization cones in 
active galaxies. 

The exponents a and j3 are variable model parameters. 
They determine the geometry and radial density distribu- 
tion. We consider a pillbox disk (a = 0), a conical disk 
(a = 1) or a flared disk (a — 9/7), see also Fig. 2] The flar- 
ing exponent o f 9/7 is chosen following; m odels for young 
stellar objects l)Chianer fc Goldreicblll997h . since there are 
no better constraints for disks around active galactic nu- 
clei. (3 determines the surface density distribution, i.e. the 
mass concentration with radius. We constructed models 
with a flat mass distribution (/3 = 0), which implies that 
the mass is concentrated in the outer regions for the con- 
ical and flaring disk. In addition, we have models where 
the mass is concentrated in the inner regions, either sim- 
ilar amounts in logarithmic intervals ((3 — —2), or less 
centrally concentrated (/? = —1). 

The radial optical depth of the disk at the equator is 
defined at 20 /im (T20), and four values are used: 5, 10, 20 
and 50. A value below 1 at 20 /im is in consistent with ob- 
served anisotropies in the mid-inf rared l|van Bemmel et all 
ll99Sll200d[Heckman et alJll992|) . The upper value is lim- 
ited by the numerical possibilities of radical. So is cho- 
sen such that the model has the desired radial equatorial 
optical depth at 20 /im. 
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Fig. 1. Left panel: the three different geometries of the dust disk. Dotted line is the pillbox model (a = 0), solid line 
the conical disk (a = 1), and the dashed line is the flared disk (a = 9/7). Right panel: an illustration of the grid setup 
in Cartesian coordinates. The grid is mirrored in the r-axis and rotated on the z-axis to achieve full spatial coverage 
of the dust disk. Refinement of the grid cells at small radii is not shown. 



The inner radius (tq) is either 3 or 30 pc, the outer 
radius (r out ) 30 or 300 pc, which results in three different 
sizes for the disk. These will be referred to as large disk 
(L = 3 — 300 pc), inner disk (1 = 3 — 30 pc) and outer disk 
(O= 30- 300pc). 

2.2. Dust characteristics 

Using the extinction efficiency tables from lLaor fc Draind 
( 1993J) we have constructed an input opacity table. The 
mass absorption coefficient n\ is obtained by summing 
the values of K\(a) for all separate grain sizes and dust 
species, using a normalized size distribution Si(a) for each 
dust species. The normalized size distribution is given by: 

a q m d ,i Xj 

Ea S *( a ) 

where a is the grain radius, md,i the mass of a dust grain 
of dust species i, Xi the fraction of dust species i and q a 
dimcnsionless exponent that can vary between —0.5 and 
—3.5. In Fig. |21 the effect of varying q on the resulting 
opacity is shown. We use only graphite and silicate dust 
in a 50/50% ratio (Xi = 0.5). The grain sizes range from 
0.001 to 10 ^m. 

For q = —3.5 the Standa rd Galactic Dust or MRN dust 
is obtained, as described bv lMathis et al. We have 

made models with q = —3.5 and q = —2.0 to test the ef- 
fect of the grain size distribution on a conical disk (a — I) 
with a constant density (f3 — 0). No silicate 10/j,m emis- 
sion feature has yet been detected in any active galaxy, 
constituting a major problem for radiative transfer models 
of obscuring tori so far. The lack of 10 /im emission might 
indicate that dust in active galaxies is dominated by larger 



grains than Galactic dust ijLaor fc Draine] | l993l). Since 
other studies point in the same direction jMaiolino et alJ 
1200111) . we have adopted q = —2.0 as standard. 

2.3. Input spectrum 

The input spectrum is defined from 10 11 to 10 17 Hz i.e. 
from 3 mm to 3nm (400 eV), and divided in 112 equally 
large logarithmic spacings. This number was chosen to 
optimize both calculation speed and frequency resolution. 
The upper limit for the frequency lies at the point where 
the dust opacity drops precipitously and the exact value 
is of little consequence. The lower limit lies far enough in 
the millimeter to ensure complete coverage of the infrared 
output spectrum. 

For the shape of the input spectrum a power-law is 
assumed: 

L(v) = Uv'J 

This allows for two more variables: the spectral slope j and 
the total luminosity. The results do not depend strongly 
on j when it is varied between and 4; we will adopt 
j = 2. The luminosity is scaled to be 10 12 L Q , integrated 
over the highest energies (100-400 eV), which is a typical 
intrinsic X-ray luminosity for radio-loud active galaxies. 
In Fig. [5] we show the spectrum that we use as input for 

RADICAL. 

2.4. Grid setup 

The radial grid is set up using tq and r out . A fixed num- 
ber of radial points is defined between the inner and outer 
radius. The grid is refined at the inner edge to prevent 
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log v (Hz) log v (Hz) 

Fig. 2. The input spectrum used for the models. On the left the flux as it is used by RADICAL, on the right the 
matching energy spectrum, using j = 2. 
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Fig. 3. The resulting opacities for different values of q, 
the exponent of the grain size distribution. From top to 
bottom q increases from -3.5 to -0.5 with steps of 0.5. 
For the models we used q — — 2 (dash-dotted line). The 
frequency resolution of this plot is higher than for the 
models. 



flux overshoots in the innermost grid cells. For the small 
disks, we have used 40 radial gridpoints; for the large disk 
60. The angular grid does not need to extend from pole 
to pole. The dust density in the wings of the Gaussian is 
negligible at a few scale heights, and the disk can be mir- 
rored in the equator. The grid is therefore defined from the 
equator to a maximum angle that extends well beyond the 
height where the opacity reaches unity. The code assumes 
that the remaining space has the same density as the grid 
cells on the edge of the defined angular raster. We use 20 



gridpoints in angular direction for all models. The result- 
ing grid is depicted in Fig. ^ 

3. The radiative transfer code radical 

In order to compute the spectrum emitted by the irra- 
diated torus, a 2-D radiative transfer calculation is per- 
formed. This is done using radical 1 , which is an all- 
purpose 2-D radi ative transfer code for axisy mmetric cir- 
cumstellar media l|Dullemond fc Turollal200(X . In this sec- 
tion we describe how we use this code to solve the transfer 
problem of a torus, partially obscuring the active nucleus. 

The dust in the torus is assumed to be in equilibrium 
with the radiation field. Each dust grain has a fine-tuned 
temperature such that it emits as much energy per second 
as it absorbs. In the case of an optically thin torus this 
temperature is easily computed, since the flux of the nu- 
cleus is known and simply goes as 1/r 2 . In such a case, no 
radiative transfer is needed, since the temperatures of the 
grains can be computed. A simple volume integral over the 
emission is then sufficient to obtain the emitted spectrum. 

Observations of radio-loud active galaxies suggest that 
the opacity of the obscuring disk exceeds unity at least up 
to mid- infrared wavelengths (~40/xm). The calculation of 
the emission changes significantly in optically thick mod- 
els. The nuclear radiation is now stopped when it is inter- 
cepted by the torus at the inner edge, and slowly diffuses 
into the dust. The inner edge is therefore expected to be 
very hot, while the interior of the torus is much cooler. The 
spectral energy distribution of such a torus with a large 
range in temperatures will be broader than expected from 



For more information, see also 

http://www.mpa-garcliing.mpg.de/PUBLICATIONS/DATA/radtrans/radicE 
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Fig. 4. Left panel: reference model compared with a model that has lower frequency resolution (dashed line). Right 
panel: comparison of the reference model (ho = 3pc, solid line) with a model that has a lower value for the inner scale 
height (ho — 1 pc, dashed line). 



an optically thin torus. Solving the temperature of the 
dust grains is much more difficult than in the optically 
thin case. This temperature depends not only on the inci- 
dent nuclear flux, but also on the infrared radiation of the 
torus itself. A grain located many optical depths deep into 
the torus will see little of the direct nuclear radiation, but 
instead will see the infrared emission of its neighboring 
grains. The radiative equilibrium of this grain depends 
therefore on the emission of the grains elsewhere in the 
torus. Every element of the torus is radiatively connected 
to every other part. 

3.1. The transfer algorithm 

At present, the best way to solve this problem is to per- 
form a 2-D radiative transfer calculation. The absorption 
and re-emission of radiative energy will be modeled as it 
seeps through the torus. The radiative transfer along each 
ray through the torus is given by: 



From this radiation field the mean intensity J„ is com- 
puted: 

J» = \J Iudn (2) 

where d£l is the solid angle. The dust temperature at every 
point can then be found by solving the equation: 

/>oo />oo 

/ B v (T)dv = \ J v dv (3) 
Jo Jo 

Once this (new) temperature structure is found, it is com- 
pared to the old one. If the relative difference is larger 
than 0.001, we go back to the integration of the transfer 
equation (Eq.[TJ , and keep iterating until the relative suc- 
cessive difference drops below the convergence criterion of 
0.001. We tested the convergence by starting the model 
with a cold (1 K) and a hot (1500 K) disk and comparing 
the resulting spectra. The results are indistinguishable, 
which indicates that the models are fully converged. 



— ^ = pK v (B v (T) — I v ) (1) 
ds 

where B U (T) is the Planck function at the local tempera- 
ture, l v the intensity, and s the path length along the ray. 
The solution of the radiative transfer proceeds as follows. 
We start with an initial guess for the dust temperature 
distribution. Using a radiative transfer integrator based 
on the algorithm of Exten ded Short Characteristics (see 
iDullemond k Turollal200(i for details), the radiation field 
at every location in the disk is computed as a function of 
frequency v and direction /i, <j> of the radiation. 



3.2. Computing the spectrum 

The spectrum is obtained from the density and tempera- 
ture distribution using a ray-tracing algorithm. For each 
discrete frequency an image of the torus is produced. 
These images are then integrated to obtain the spectrum 
F v . To make sure that the pixels of these images accu- 
rately sample all relevant scales in the system, we dis- 
tribute the pixels of the image not in a rectangle (as in 
normal digital images), but in concentric circles around 
the central source. Each circle contains 64 pixels, and the 
circles themselves have radii such than for every radial 
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og A (/im) log A (/^m) 

Fig. 5. The reference model, using q = —2, a = 1, j3 = 0, T20 = 10, with a disk size from 3 to 30 pc. Different lines 
are for different aspect angles, the top line corresponding to 9 = 5°, next are 30, 50, 60, 70 and 90°. In the right panel 
a zoom on the infrared regime is displayed. 



gridpoint on the computational domain there are two cir- 
cles in the image. In this way there will be an equal number 
of pixel-circles per logarithmic interval in radius, and this 
guarantees that both the very inner parts of the disk (few 
pc) as well as the far outer parts (tens to hundreds of pc) 
will be sampled properly. 

The spectra can be computed at any requested incli- 
nation. We select inclinations (9) of 5, 30, 50, 60, 70 and 
90 degrees. For small 9 the inclination cannot be smaller 
than 0.1 radian, due to limitations of the code, so its low- 
est value is 5.7 degrees. All 2-D/3-D effects are included. 
For edge-on inclination (9 = 90°) one will see only the 
outer edge of the disk, while at more face-on inclination 
(say 9 — 5°) one will also see the central source and the 
warm inner regions of the disk. The angle at which the in- 
ner regions are shielded from view depends on the opening 
angle of the torus and the optical depth in our models. 

3.3. Flux conservation and resolution tests 

In the ideal situation, the torus will emit exactly the same 
amount of luminosity as it has absorbed from the central 
source. This requirement of flux conservation is a pow- 
erful check on the results of the code, and this check is 
in fact always performed by radical at the end of each 
run. Moreover, radical is built in such a way that radi- 
ation in optically thin areas of the system does not loose 
or gain flux through numerical diffusion. Unfortunately, 
in areas of large optical depth, flux conservation is not 
intrinsically guaranteed by the code, and occasionally the 
flux can be slightly over-predicted. In most cases we have 
energy conservation to within about 10%, but sometimes 
(larger disk and/or high opacities) this error can go up to 
20%. Refinement of the grid and fine-tuning of the bound- 
aries of the computational domain usually help to reduce 
this error, but there are feasibility limits to the number of 
grid-points that can be used. In this paper all models stay 



within this 20% tolerance, and most of the models in fact 
stay within ~10%. 

The frequency resolution has to be high enough to pre- 
vent flux losses or overshoots. In Fig.^Jthe effect of a lower 
frequency resolution is shown: with half the number of fre- 
quency points the flux deviates less than a factor of two. A 
higher frequency resolution does not significantly change 
the results. 

Changing the scale height at the inner radius also af- 
fects the output spectrum. As shown in Fig. 0] decreasing 
ho mainly results in a lower output flux in the far-infrared. 
Since the disk is thinner, there is less cool dust in the disk. 
The decrease of flux is proportional to the decrease in ho . 
It only marginally affects the shape of the spectrum, there- 
fore it can be used to tune the results to observations for 
individual objects. When comparing colours, the assumed 
value of ho is irrelevant. 

3.4. Improvement on previous methods 

The large parameter space probed by our models is unique 
and allows for a full comparison of observations with these 
models. The dependence of the output spectrum on disk 
parameters has not been studied in this much detail be- 
fore. For simplicity we have chosen to deal with a smooth 
dust distribution, before assessing the effects of a clumpy 
medium. One of the main problems has always been to re- 
produce the broad infrared SED observed in active galax- 
ies: the effects of the geometry of the torus, its density, 
size, or opacity are not clear. We compare the models to 
available broad band data in order to assess the width of 
the resulting spectrum. Instead of fitting a large sample of 
objects, we will compare the colours of our models to the 
observed colours. Only for Cygnus A, a powerful nearby 
radio galaxy, will we study the actual fit in detail. 

A second important issue has been the occurence of 
10 u,m emission at small inclinations in most literature 
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Fig. 6. The temperature and density distribution for the conical disk (left two panels) and the pillbox disk (right two 
panels). Upper panels are in z coordinates, lower panels in 9/tt. Grey scale shows the density distribution of the dust, 
the solid lines are temperature contours, corresponding to 600 K (leftmost) and decreasing with steps of 50 K. 



models, which is not observed. Although our frequency 
resolution does not allow an in depth study of the be- 
haviour of the 10 /im feature, we will discuss the effects 
of altering the grain size distribution on our models. The 
powerful radio-loud active galaxies lack proper observa- 
tions of the 10 fim region, and since we mainly aim to un- 
derstand these objects, no direct comparison to infrared 
data will be presented in this paper. In a consecutive pa- 
per, high resolution models will be discussed, allowing 
for a direct comparison with new observations from the 
SIRTF satellite for a range of active galaxy types. 



4. Results 

In order to understand the effect of the different variables, 
we have defined a reference model, to which all other mod- 
els are compared. The relevant parameters are given in 
Tabled The reference model is chosen to be a conical disk 
with constant surface density, a size from 3 to 30 pc and 
optical depth of 10. This will allow a direct comparison 
to literature models, which generally have similar charac- 
teristics. However, we will use q = —2.0, which docs in- 
troduce a different behaviour when considering the 10 /im 
silicate feature, but which is a more realistic value when 
considering active galaxies. The optical depth of 10 is cho- 
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Fig. 7. Comparison of the models with different varying parameters for a fixed inclination of 9 = 30°. The non- varying 
parameters have the values as for the reference model: a — 1, (3 — 0, T20 = 10 and size 3-30 pc. 



sen to account for observed differences at 40 /xm restframe 
wave l ength in quasars and radio gal axies ijHeckman et al.| 
Il992tlvan Bemmel et al1ll998ll2000|) . The size of the disk 
corresponds to sizes found in Seyfert galaxies, although it 
may be significantly larger in powerful radio galaxies. 

4.1. The reference model 

The model spectrum is shown in Fig. The spectrum 
peaks between 10 and 20 /xm, depending on inclination. 
For the different inclinations there are clear flux differ- 
ences with inclinations (anisotropies) up to 100 /xm. The 
knee at 110 /xm is a dust opacity feature. There is a slight 
indication of 10 /xm silicate absorption at larger angles. 

The temperature and density distribution for the ref- 
erence model are depicted in Fig. El The temperature in 



the disk is lower than in its optically thin surroundings, 
which is to be expected. In the inner regions the temper- 
ature drops steeply with radius, flattening at larger radii. 
Within the disk, the temperature in angular direction is 
more or less constant, with a steep increase at the edges 
to reach the temperature of the surrounding medium. 

Below is described what the effect of varying the differ- 
ent parameters has on the model results, which is shown in 
Figs. IHland El These figures are a concise representation 
of the effects of changing the parameters. A complete set 
of hundreds of figures has been used to study the effect of 
the variables on anisotropy and 10 /xm absorption in the 
models. 

The depth of the 10 /xm absorption sh own in Fig. 
is det ermined using a similar method as Pier & Krolik 
a power-law is fitted through the 6.5 and 13.9 xtm 
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Fig. 8. As Fig.0 but for inclination angle = 70°. 



points. However, since the model curves have a grey-body 
shape, fitting a power-law will mimick silicate emission 
where none is present. This is circumvented by adding 
to the power-law fits the residual from the model with 
5° inclination, where no silicate emission or absorption 
occurs, resulting in a continuum that approaches a grey- 
body shape in the region 5-15 /urn. Subsequently, the dif- 
ference between the continuum and the 9.7 /im gridpoint 
is measured. The depth rg.7^ m is positive for aborption. 

— a: changing the physical shape of the disk to a pillbox 
(a = 0) or flaring disk (a = 9/7). The flaring disk and 
conical disk have similar spectra, slightly varying with 
inclination. The pillbox disk is less luminous at all in- 
clinations, but this is more obvious at smaller angles. 
It is hotter than the reference model, with the peak at 
shorter wavelengths for smaller inclinations. The wave- 



length of the emission peak does not shift significantly 
with inclination, but there is an indication of increased 
anisotropy. The flared models have a higher degree of 
isotropy than the reference model, but show a simi- 
lar shift in the peak wavelength with inclination. They 
also show a hint of 10 /urn absorption at small inclina- 
tions, but this effect depends strongly on the optical 
depth and the size of the disk. 
— (3: changing the surface density from a constant to a 
more centrally concentrated configuration (f3 = —2). 
For fixed a, the models differ slightly at large inclina- 
tions, and marginally at small inclinations. The more 
centrally concentrated, the more the spectrum moves 
to shorter wavelengths. At large 8, there is a hint of 
10 /im absorption, which deepens when (3 changes from 
to -2. The depth does not depend on the value of T2o- 
The 10/Ltm absorption is absent in the 30-300 pc disk. 
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Fig. 9. Comparison of the 10 /urn silicate absorption in the reference model when using different opacity laws. The 
solid line represents MRN dust (q = —3.5, Galactic), the dashed line the larger grain dust (q = —2) adopted in this 
paper. Left panel: A zoom on the spectrum for inclinations 5°, 30°, 50°, 60°, 70° and 90° (from top to bottom). Right 
panel: Depth of the 10 /im absorption feature with inclination (positive values indicate absorption) . 



— T20: changing the optical depth at 20 /im from 5 to 
50. The higher opacity models have broader infrared 
spectra and shorter peak wavelengths at small angles. 
For larger inclination the high opacity models are not 
much broader, but clearly have more far- infrared emis- 
sion. The peak wavelength hardly changes at larger an- 
gles. An important effect is the increase of anisotropy 
with increasing T20, but the wavelength at which the 
dust becomes optically thin does not vary significantly. 
There is a hint of 10 /im absorption in the low opacity 
models at large inclination. 

— size: the outer disk emits at much longer wavelengths 
than the inner disk, and shows a stronger anisotropy at 
longer wavelengths. The large disk has a much broader 
spectrum at small angles, indicating a larger temper- 
ature range. At large inclinations the outer and the 
large disk are almost identical. The large disk has clear 
10 /im absorption. A hint of this is seen in the inner 
and outer disk. 

— q: using Galactic dust instead of a large-grain dom- 
inated mixture. For larger inclinations the 10 /im ab- 
sorption becomes deeper. There is a hint of 10 /im emis- 
sion, which strongly depends on the optical depth and 
disk size. For smaller optical depth, the emission is 
clearly present. For larger optical depth, the absorp- 
tion is only seen in intermediate inclinations. A direct 
comparison of 79. 7^ shows minor silicate emission at 
inclinations around 30° , but up to a factor of 5 increase 
of silicate absorption at angles larger than 40°. In the 
outer disk, the absorption is deeper at larger inclina- 
tions, but there is a clear emission feature for small 
angles. The anisotropy is less pronounced than in the 
reference model, although it is visible up to 100 /im as 
well. The peak wavelength is unchanged. 



It is illuminating to examine the spatial distribution of 
the emission in the disk. In Fig. ^| we show images of the 
reference model, at 30 and 200 /im, and with inclination 
of 30 and 70°. The images display the entire angular grid, 
causing them to have a circular cut-off, independent of 
the inclination of the disk. At 30 /im the disk is optically 
thick and only the inner regions radiate. At 70° the inner 
regions are not seen directly, but some emission is seen 
from the opposite inner side of the torus. At 200 /im the 
dust is optically thin and the entire disk will radiate. At 
30° the outermost regions are not as bright as the rest, 
since they are somewhat cooler. At 70° the inner regions 
are now seen direct, since the dust is transparent. They 
are the brightest part of the disk, since this is where the 
column density is highest. 

The effect of varying the parameters is shown in 
Fig. 1111 The top row shows again the reference model. 
The second row shows the pillbox models (a = 0), which 
are clearly less bright. This difference is most pronounced 
at 200 /im. They are of course flatter at 70° . The third row 
displays the (3 = —2 models, i.e. models with equal mass 
per logarithmic interval. Here the inner regions dominate 
all the images, since they represent the densest parts of 
the disk. This is especially clear in the 200 /im images, 
where the disk is optically thin. The final row contains 
the images for a disk with T20 = 50. At 30 /im the emis- 
sion is more centrally concentrated and the contrast with 
the outer regions is larger than for the reference model. 
However, for the 200 /im images the opacity is irrelevant, 
and the images are comparable. 

4.2. Comparison with models in literature 

The resulting spectrum from our reference model can be 
compared to literature models with comparable disk and 
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Fig. 10. Images of the reference model, using a circular aperture. Top panels are for 30° inclination, bottom panels 
for 70°. Left panels are the 30 /jm images, right panels are at 200 /im. The grey scale is inverted: dark areas represent 
more emission. The central white region is where the disk ends and the central emission source is located. The outer 
circle represents the aperture. 



dust parameters. The main difference between our model 
and the literature models is the absence of significant 
10 /im emission at small inclinations, and only weak 10 /im 
absorption at larger inclinations. However, since all litera- 
ture models use Standard Galactic Dust, a comparison at 
10 /im is unrealistic. 

A qualitative comparison of our reference model with 
literature models with comparable parameters reveils 
that the general shape of our outpu t spectrum is best 
matched by the iPier fe Krolikl lll992l) models (he reafter 
PK92). The models bv iGranato fc Danesd l|l994|) (here- 
after GD94) also match the shape q uite well, but the 
lEfstathiou fc Rowan-Robinsonl l)l995|) (hereafter ER95) 
models are very different, with a much flatter SED in 
the mid-infrared and often a peak in the near-infrared. 
Our reference model becomes optically thin longward of 
100 /im. Only the PK92 models match this, other mod- 
els all become optically thin between 60 and 80 /im. In all 
literature models changing the inclination causes a con- 



tinous decrease in flux at a fixed wavelength. The GD94 
models show a very discrete change in flux between 45 and 
75 degrees inlination. As a consequence, there is no obvi- 
ous difference between their models with and 45 degrees 
inclination, and neither between the models with inclina- 
tions of 75 and 90 degrees, except for the very thick and 
compact geometries. 

For a more quatitative comparison we have calculated 
an anisotropy value, defined by taking the flux ratio for 
face-on over edge-on oriented disks at 20 /im. For the liter- 
ature models that have parameters closest to our reference 
model, we find this ratio to be ~ 5. This is consistent with 
the anisotropy value in our reference model. 

Cl umpy torus models, as described bv iNenkova et all 
(2002) are in general flatter than our reference model, 
which is especially prominent at small inclinations. Their 
anisotropy value is also ~ 5, as in our reference model, and 
they do not become optically thin shortward of 100 /im. 
Depending on the mean free path of the clumps, the mod- 
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Fig. 11. Images of the disk when varying the different parameters. The top row is the reference model, then from 
top to bottom per row are changed: a = 0, j3 = —2 and T20 = 50. From left to right the images are for 30° at 30 and 
200 /im, and 70° at 30 and 200 /im respectively. 



els have a more pronounced peak around 60 /im, and be- 
come more luminous. 

A more detailed comparison is not feasible, since the 
models in literature use different setups and calculation 
methods. All observed differences between our results and 
literature models can be ascribed to this. 



5. Understanding the behaviour 

The output of a dusty disk depends mainly on two proper- 
ties: the temperature and the mass of the dust. If the dust 
is further away from the emission source, it will be cooler, 
and therefore emit at longer wavelengths. At shorter wave- 
lengths the dust is self-shielding, since the dust is optically 
thick for its own radiation. The inner regions will become 
very hot, but these can only get rid of the energy at wave- 
lengths where the dust is optically thin. As a consequence, 
the inner regions will be the brightest at the shorter wave- 
lengths, and only in the far-infrared will the entire disk 
start to radiate. This effect is clearly seen in the images 

(Fig.nuj. 



In a pillbox model, the r = 1 surface is 'squeezed' 
more to the equatorial plane and, hence, compared to the 
other models less radiation of the central source is inter- 
cepted. Pillbox models are therefore less luminous. The 
difference between conical and flaring disks is less pro- 
nounced. Nevertheless, for a flaring disk, some dust at 
large radii will have a direct line of sight to the central 
source and be warmer. This effect is best observed at large 
inclinations (see Fig- EJ - 

When we vary (3 from to -2, more mass will be con- 
centrated in the inner regions. This dust will be hotter, 
and therefore brighter in the near- and mid-infrared. Since 
the dust is more centrally concentrated for (3 — —2, there 
is less mass in the outer regions, and therefore less cool 
dust, as seen in Fig. ^2 As a consequence the far-infrared 
emission will decrease in flux, and the entire spectrum 
shifts to shorter wavelengths. 

Increasing the optical depth increases the amount of 
absorbed radiation. For higher values of T20 it will be more 
difficult to transport absorbed energy into the torus. As 
a consequence the inner parts get hotter, but the outer 
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Fig. 12. Colour-colour plots for models with /3 = 0, T20 = 10 and varying a and disk size. Lines indicate varying 
inclination: from left to right corresponds to edge-on to face-on orientations. Solid lines are for the inner disk, dotted 
for outer disk, dashed for large disk. The objects are coded as follows: open circles=quasars, open squares=narrow-line 
radio galaxies, open triangles=broad-line radio galaxies, filled circle=Heckman quasar sample, filled square=Heckman 
radio galaxy sample, filled triangle=Cygnus A, crosses=Seyfert 1, plusses=Seyfert 2, asterisks=ultra- luminous galaxies, 
open stars=star burst galaxies. 



parts will remain cooler. This results in a wider tempera- 
ture distribution, which translates into a broader energy 
distribution for small inclinations (Fig. 0). For large in- 
clinations, the increased amount of cold dust is most im- 
portant. Since the inner regions are shielded from direct 
view, the energy distribution shifts to longer wavelengths. 

The distance to the central source is critical in the 
models. At large inclinations, we cannot view the inner 
edge of the torus direct. The hot dust is shielded from view 
by the cool outer regions. For this reason, the large disk 
and outer disk have comparable spectra. Overall the inner 
disk is much warmer and brighter at the shortest wave- 



lengths. At small inclination, the whole disk is visible and 
both the inner and large disk, which have hot dust close 
to the central source, are much brighter at shorter wave- 
lengths. The large disk and the outer disk show very sim- 
ilar behaviour at the longest wavelengths where the dust 
is optically thin. The large disk produces less near- and 
mid-infrared emission than the inner disk, even though its 
inner radius is similar, since the temperature gradient is 
steeper. This is also related to the model setup: since the 
disk is larger, a smaller value for £0 is required to obtain 
the requested equatorial optical depth, and hence the total 
opacity is lower. 
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Fig. 13. Colour-colour plots for models with a = 1, T20 = 10 and varying /3 and the disk size. Lines and symbols as 
in Fig. U3 



Assuming q = —2 instead of using Galactic dust 
mainly affects the 10 /jm silicate feature. Although our 
models clearly lack significant silicate emission, they still 
display clear absorption, albeit not as deep as the Galactic 
dust models. The detailed behaviour of the silicate feature 
will be the purpose of future high frequency resolution 
models. 

6. Comparing the models to observations 

6.1. Colour comparison 

Using the known infrared colours of samples of active 
galaxies, we can attempt to constrain the parameters of 
the obscuring medium. We will use the IRAS bands (12, 
25, 60 and 100 ^m) which are also available as ISO bands. 



We have select ed all 3C radio-loud a c tive galaxies from 
the l iterature l|Golombek et al.1 Il988t iNeugebauer et all 
1986) with a least three solid detections in three of the 
IRAS bands, and with small redshifts (z < 0.5). It should 
be mentioned, however, that IRAS detected only a few 
radio-loud active galaxies, which all lie close to the de- 
tection limits of the instrument. In addition we selected 
a sample of starburst, ultra-luminous and Sey fert galax- 
ies with four IRAS de tections from literature l)Sekiguchi 
ll987tlMa,rKeTrtvlll98^ . 

For all radio-loud objects we as sured that the infrared 
emission is of thermal nature. From lHeckman et alJ (^.992) 
we have constructed the colours for samples of quasars 
and radio galaxies. All data are overplotted on the mod- 
els, while varying different parameters in Fig. 1121 1131 and 
n~T1 The a symbol on the axis, annotated with two wave- 
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Fig. 14. Colour-colour plots for models with a = 1, p = and varying T20 and disk size. Lines and symbols as in 
Fig. El 



lengths, represents the spectral index or colour, defined as 
F oc v a , while the a symbol in the plots, without anno- 
tation, refers to the geometry of the disk. We will refer 
to the colours as hot (a^|), warm (af"), and cold (a^ ). 
For an increase in spectral index the flux at the shorter 



wavelength increases with respect to the flux at the longer 
wavelength. Hence an increase in the value for the colour 
implies more flux at shorter wavelength. 

The objects within a given class cluster together. The 
quasars have a narrow range of hot and cold colours, both 
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lie between and —1. The warm colour ranges from 0.3 to 
—2. The narrow- line radio galaxies have cooler hot colours 
than the quasars, ranging from —0.5 to —2. The warm 
colours show approximately the same values, from to 
— 1.5, but the cold colour has a large range from 1 to 
almost —2. For the broad- line radio galaxies the hot colour 
is comparable to that of the quasars, but the warm colour 
is clearly warmer, having typical values between 1 and 
—0.5. The cold colour is only determined for one object, 
the rest lack 100 /xm detections. The starburst and ultra- 
luminous galaxies lie clearly in a different region than the 
active galaxies in the hot vs. warm plot, while their cold 
colours are comparable to the cold colours of the radio- 
loud active galaxies. The Seyfert 1 colour distribution is 
indistinguishable from the colours for the radio-loud active 
galaxies, but the Seyfert 2 galaxies overlap both the active 
galaxies and the starbursts in their hot and warm colours. 

There are no objects which have a positive hot colour 
in the entire sample, and only a few which have a pos- 
itive warm colour. Nevertheless, the colours of the three 
radio-loud classes do not separate out too much. From 
an infrared point of view, these three classes are a some- 
what arbitrary classification of a continuum distribution 
of types. 

6.1.1. Models with (3 = and r 20 = 10 

In Fig. 1121 we have fixed (3 and T20 , and varied the disk 
size (3 different lines in plot), a (the different plots) and 
the orientation (along the lines, edge-on to face-on is from 
left to right). The overplotted data are for objects which 
are presumed to represent a continuum of different ori- 
entations, from quasars to narrow-line radio galaxies. The 
filled symbols are for ensembles of quasars and radio galax- 
ies. 

Changing a from a — (pillbox) to a = 1 (conical) has 
a profound effect on the colours of the disk. A further in- 
crease to a = 2 (flaring disk) has little effect. These effects 
are in line with the SEDs shown in Figs. and |H1 

6.1.2. Models with a = 1 and T20 = 10 

In Fig. 1131 the value for a is fixed, as is T20 = 10. 
Decreasing (3 to — 2 the models move up, except the outer 
disk, indicating more positive warm and cold spectral in- 
dices. The hot colour only increases with decreasing (3 in 
the large disk. 

6.1.3. Models with a = 1 and (3 = 

So far we have fixed the value of T20, but this can of course 
also vary. In Fig. El we have fixed a = 1 and (3 = and 
varied T20 in the different plots. Along the lines we still 
vary the inclination. 

Increasing T20 will widen the range of the hot and warm 
colours of the models with changing the inclination. The 
warm colour also becomes cooler with increasing opacity. 



The cold colours have approximately the same range when 
varying T20, but become cooler for higher values. 



6.1.4. Comparison with observations 

As judged from the plots, for all models the cold colours 
appear inconsistent with the observations. This implies 
that either the observed 60 /im flux or the 100 /xm flux, or 
both arc not fitted by the models. The hot colours in the 
models span the same range as the observed colours for 
active galaxies, but in the hot vs warm colour plot, the 
warm colours in the models do not cover the entire range 
of observed values. This indicates that already at 60 /im 
the spectrum starts to deviate. Hence, both the 60 and 
100 /im emission in the active galaxies are dominated by 
a different component. 

At shorter wavelengths, the difference in hot and 
warm colours between starbursts and active galaxies in- 
dicates that the 12 and 25 /im emission in active galaxies 
arises from a component unrelated to star-formation, 
most likely the active nucleus. Longward of 60 /im the 
colours of active galaxies are indistinguishable from the 
colours of non-active galaxies, which implies that the 
emission process is the same in both groups. This favours 
star-formation as the dominant process of far-infrared 
emission in active galaxies. A similar conclusion is 
reached by studies of infr ared colours o f Seyfe rt galaxies 
longward of 100 /im by ISpinoglio et al. (2002), who find 
that this emission is indepe ndent of the prese nce of an 
active nucleus. Furthermore. IPrieto et alJ I 200lh conclude 
that coronal line strength is strongly correlated with 
near- and mid-infrared emission, but not at all with 
far-infrared, indicative of a different heating mechanism 
for the far-infrared emission. 



In the hot vs warm colour plot the narrow-line and 
broad- line radio galaxies are well fitted by the models. 
For the other colour plots only a few narrow-line ra- 
dio galaxies are still within the range of model colours. 
The outlying quasar in the hot vs warm colours is 3C 48, 
known to have a large star-formation rate, and possi- 
ble all its infrared emission arises in star-forming regions 
l)Canalizo fc Stocktonll200o|k Its warm colour is compara- 
ble to the warm colour of objects with known large star- 
formation rates, but its hot colour is clearly warmer, in- 
dicative of very hot dust. This indicates that a significant 
contribution from an active nucleus will cause objects to 
have a more positive warm colour, while their cold colours 
are indistinguishable from the colours of star-forming ob- 
jects. 

The inner disk does not fit the observed colours well, 
except in the hot vs warm colour plot. The large and 
outer disk do better for a = 1 and a — 2 disks and 
large inclinations, but also fail in the hot vs cold colours. 
The pillbox disks do not fit any of the observations. For 
decreasing f3 all models move away from the observations. 
Increasing t 2 q moves the models closer to the observed 
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colours, but the discrepancy in cold colours is not 
explained. 

We note with interest that the ultra-luminous galax- 
ies are well fitted by the high opacity models in the cold 
vs. hot and cold vs. warm colour plots. Also Cygnus A 
lies closer to the models than the remainder of the ra- 
dio galaxy population. For the ultra-luminous galaxies it 
is known that star-formation dominates the far-infrared 
output, and also Cygnus A is has strong evidence for 
significant star-formation in the host galaxy. Hence, star- 
formation can mimic the SED of an obscuring torus and 
its relative contribution must be assessed by other means, 
e.g. PAH emission, optical absorption lines, to fully disen- 
tangle the two infrared components 

6.1.5. Conclusions from colour-colour plots 

Most objects have negative cold spectral indices which 
are not explained by the models. This could very 
well be due to a contribution from other components, 
such as additional cold dust from star-formation re- 
gions. Observational evidence for a second dust com- 
ponent heated by stars in active galaxies is found in 
several studies llPrieto et alJ 200ll ISninoriio et~al]l2002l 
van Bemmel fc Bar thcl 2001|). Therefore, we conclude 
that the torus emission only dominates the SED to about 
25 /im, while the observed 60 and 100 /im flux arise in a 
second infrared component. 

With respect to the unification theory, there is no 
trend in the data for radio galaxies and Seyfert 2 to be 
more consistent with large inclination models, or quasars 
and Seyfert 1 with small inclinations. In the hot vs warm 
plots all radio galaxies lie closer to the models than the 
quasars, while the objects with known large star-formation 
rates occupy a clearly distinct region, having cooler hot 
and warm colours. The broad-line radio galaxy colours 
are consiste nt with the lack of large scal e cold dust, as 
discussed in Ivan Bemmel fc Barthe If the torus is 

the only infrared component in broad-line radio galaxies, 
they should also be consistent with the models in their 
cold colours, but we lack detections at 100 (im to confirm 
this. As emphasized above, the differences between mod- 
els and observations are indicative of a second component 
in the far-infrared. 



6.2. Comparison with Cygnus A 

For individual objects, the model parameters can be con- 
strained by a direct fit to the observed infrared spec- 
trum. One of the best objects to test our models is the 
narrow-line radio galaxy Cygnus A, since it is the best ob- 
served radio-loud active galaxy in the infrared, and has a 
rather well determined central continuum. IRAS and ISO 
have observed Cygnus A in the infrared and additional 
(sub-)millimcter data are available to determine the non- 
thermal contribution. From optical observations the size 



of the obscuring torus has been constrained to lie between 
a few parsecs inner radius, up to maximum of a fe w hun- 
dred parsecs outer radius (jlmanishi fc Uenoll2000h . Also, 
we know that the angle of orientatio n of the radio jet is 
~ 50 degrees (|Carilli fc Barthelll996h . allowing us to vary 
9 around this value. For the models we fix a = 1 and 
,9 = 0, since these are the best values from the colour- 
colour fitting. We use q = — 2 for the dust grain size dis- 
tribution. 

We have constructed models normalized to the 
observed hard X-ray luminosity and spectral slope, 
and taken into account size limitations of the torus in 
Cygnus A. The inner and outer radius are varied with 
one extra value: r$ — 3, 10 and 30 pc, r out — 30, 100 and 
300 pc. Models are constructed for the same four values 
of T20 and spectra are obtained for varying inclinations. 

In Fig. we show several fits to the spectrum of 
Cygnus A, with different variables. A small scaling factor 
is needed to fit the models, which is either due to a dif- 
ferent intrinsic luminosity, or different value of ho, which 
are both kept fixed. 

6.2.1. Tuning the fit: size of the disk 

The size of the disk has the largest effect on the width 
of the SED. Good fits are obtained for all optical depths 
for a disk with r = 10 — 100 pc. The reference model does 
not fit the SED, as expected from the colour-colour anal- 
ysis. It overproduces near-infrared emission, and cannot 
account for all the far-infrared emission. The large and 
outer disk can fit the far-infrared output of Cygnus A, but 
both underpredict the observed mid-infrared flux. The in- 
ner radius of the disk can be carefully tuned to the shape 
and intensity of the mid-infrared observations. For an in- 
ner radius of 10 pc the mid-infrared is well fitted by the 
model. 

With an outer radius of 30 pc all four mid- infrared 
points are well fitted by the different opacity models. The 
model with an outer radius of 100 pc fits all the data, 
although the fit in the mid-infrared is somewhat worse. 
The longest wavelength far-infrared point is only fitted by 
the 10-100 pc disk with T20 = 5. Increasing the outer ra- 
dius to 300 pc improves the fit for other opacities at the 
longest wavelengths, but then the fit in the mid-infrared 
deteriorates. The models with an outer radius of 100 pc or 
more peak at 60 /im, as observed, and the slope of the far- 
infrared and sub-millimeter tail is reasonably consistent 
with the observations. 

We have no problem reproducing the width of the ob- 
served SED with a single large torus model, as did previ- 
ous models. However, from optical imaging and spectropo- 
larimetry we kn ow that Cygnus A has large scale dust an d 
star-formation IjJackson et al.l l 1998t iFosburv et aflll999fr . 
which is expected to contribute to the far-infrared emis- 
sion. From the colour-colour analysis we also conclude that 
the far-infrared emission beyond 60 /im is generally not 




Fig. 15. Upper left panel: model fits of a conical disk to the observed SED of Cygnus A. The error bars on the observed 
luminosities are smaller than the plot symbols. The model disk extends from 10 to 100 pc from the nucleus and has an 
inclination of 50°. Lines are for different opacity models (T20): 5=dotted, 10=dashcd, 20=solid and 50=dash-dotted. 
Upper right panel: as left panel, but for a disk from 10-30 pc. Lower left panel: fits of disks with T20 = 10 and inclination 
of 50 degrees. Lines arc for different sizes: dotted=10-100pc, dashed=10-30pc, solid=3-30pc, dash-dot=30-300pc, 
dash-triple dot=3-300pc. Lower right panel: fits of disks with T20 = 10, and size 10 to 100 pc. Different lines are for 
different inclinations. From top to bottom: 30 (dotted), 50 (dashed), 60 (solid) and 70 (dash-dotted). All panels have 
a = 1 and (3 = models. 



produced by a torus. Therefore, the 10-30 pc disk provides 
the best model for the torus in Cygnus A. 

6.2.2. Tuning inclination and opacity 

Small variations in inclination have little effect on the re- 
sulting SED and models with inclinations between 5 and 
60 degrees fit equally well. Models with higher inclinations 
are ruled out. 

The opacity cannot be fully constrained from the mod- 
els. For the 10-100 pc disk the t 2 o = 10 and 20 models fit 
best, which is also the case in the 10-30 pc disk. However, 
in the latter the opacity is not constrained by the far- 



infrared, which leaves only four datapoints to optimize 
the fit. The effect of varying a and (3 has been described 
above, and is similar in the Cygnus A models. Changing 
them will not significantly improve the fit. 

6.2.3. Conclusions from the Cygnus A fit 

From the observations of Cygnus A we can constrain the 
size and inclination of the obscuring torus. The geometry 
and density distribution are fixed, varying these parame- 
ters would not significantly change the output spectrum. 
Assuming significant star-formation in the host galaxy 
produces the far-infrared emission, the best fit for the 
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torus is obtained for a disk with size of 10-30 pc, incli- 
nation of 50 degrees and optical depth of 10-20. A single 
large disk fits all infrared points, but does not allow a 
contribution from the observed large scale dust, and is 
therefore not fully consistent with optical observations of 
young stars and large scale dust in Cygnus A. 

7. Conclusions 

Radiative transfer models have been constructed to study 
the emission of a dusty torus around an active galactic nu- 
cleus. We have varied optical depth, disk geometry, disk 
size, density distribution and dust grain size distribution, 
to study the effect on the resulting SED. In most models, 
the dust becomes optically thin around 100 /zm, which im- 
plies that the emission is isotropic at longer wavelengths. 
The models are compared to the observed infrared colours 
of a sample of radio-loud active galaxies, and furthermore 
fitted to observations of Cygnus A. The problem with 
previous models was often that the resulting SED was 
not broad enough. We have identified this as a matter of 
size and/or opacity: the larger the disk, or the higher the 
opacity, the broader the temperature distribution, and, 
as a consequence, the broader the spectrum. In the fit to 
Cygnus A we have no problem to reproduce the entire 
width of the observed infrared SED. 

From the comparison to observed colours, we conclude 
that the tori in active galaxies are relatively large. The 
inner disk model does not fit any of the observations: the 
inner radius more likely lies around 10 pc. The outer radius 
is not well constrained. Pillbox models are rejected on ba- 
sis of the colour-colour analysis, both the conical and flar- 
ing disks are possible. The surface density seems constant 
with radius; for more centrally concentrated models the 
colours are not consistent with observations. Varying the 
optical depth can sometimes improve the level of agree- 
ment, but its effect is diverse. An optimum value cannot 
be derived. 

The objects that are best fitted by the models are the 
radio galaxies. In the hot vs warm colours all broad-line ra- 
dio galaxies fall entirely in the range of the model colours. 
This is consistent with observational evidence that in 
other objects a second dust component contributes, which 
is not related to a central obscuring torus. This second 
component alters the shape of the spectrum, and hence 
the colours, which explains the inconsistency of the mod- 
els with the observed cold colours. We conclude that the 
torus dominates only at 12 and 25 /im. As a consequence, 
selecting active galaxies is best done using their 12 and/or 
25 /zm power. However, at these wavelengths the quasars 
are expected to be brighter, for which bias one has to cor- 
rect. 

Eel et al find a relation between the 60 fjxa lu- 

minosity and radio power of active galaxies: the more pow- 
erful active galaxies are more luminous at 60 /im. They 
interpret this as evidence that the 60 /im emission is dom- 
inated by reprocessed radiation from the obscuring torus. 
However, the above analysis suggests that the 60 /*m emis- 



sion in active galaxies is dominated by a second compo- 
nent, and hence its power is not indicative of the power of 
the active nucleus. If star-formation dominates the 60 /zm 
emission in active galaxies, the correlation would imply 
that the more powerful active galaxies have higher star- 
formation rates. A better understanding of the process 
dominating the 60 /zm emission can be obtained by com- 
paring individual fits to objects with and without star- 
formation, such as 3C48 and the broad-line radio galax- 
ies. 

For Cygnus A, the inner radius of the disk is con- 
strained to lie around 10 pc. This value is found to be 
most consistent with the observations in the colour-colour 
plots. The outer radius can be anywhere within 30-300 pc, 
but the best fit is obtained with a value of 30-100 pc. 
All values are within the observational limits set by opti- 
cal studies. Although a single torus model fits the entire 
infrared SED, there is strong observational evidence for 
a second dust component in Cygnus A. In addition, the 
quality of the fit in the mid-infrared improves somewhat 
for a smaller disk. An outer radius of 30 pc is preferred to 
allow for the observed star-formation to contribute to the 
far-infrared emission in Cygnus A. 

The consequences for unification models are ambigu- 
ous. On one hand, we can fit the entire SED of Cygnus A 
with a single torus model, while we know that other pro- 
cesses must contribute to its infrared emission. This is in 
itself evidence in favour of the presence of an obscuring 
torus. On the other hand, the cold colours of most ob- 
jects are inconsistent with the torus models, indicating 
that a second component must be present, and a good fit 
to the entire SED is not possible. We know that this sec- 
ond component must also be present in Cygnus A, and 
therefore, the relative contribution of star-formation must 
be assessed with other means to thoroughly disentangle 
the infrared SED of active galaxies. Only with a proper 
decomposition of the SED can the actual emission from 
an obscuring torus be studied in detail. 

8. Future work 

Infrared observations provide a powerful tool to study the 
emission from an obscuring torus in any type of active 
galaxy, especially at near- and mid- infrared wavelengths, 
where the emission from the torus is observed directly. The 
models presented in this paper are a next step in increas- 
ing our understanding of such tori. However, improvement 
of the models can be expected on several fronts. 

As mentioned, the 10 /im absorption feature is sensi- 
tive to the parameters of the model. Unfortunately, no 
data presently exist for radio-loud active galaxies, due to 
their faintness at mid-infrared data and lack of capable 
instruments. SIRTF is expected to have the sensitivity to 
observe a large sample of radio- loud active galaxies, over 
a large range of rcdshifts. By constructing models with 
higher frequency resolution around 10 /im, predictions can 
be made for future observations of 10 /im emission or ab- 
sorption in radio- loud active galaxies. In addition the spec- 
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tral index of the mid-infrared spectrum is sensitive to the 
parameters, especially a^, and to some extent o^s- Since 
this is the wavelength range where we are reasonably sure 
we observe emission from the torus, high resolution spec- 
tra arc needed. The SIRTF satellite will be able to obtain 
spectra up to 90 /im for active galaxies over a large range 
of redshifts, which will provide a good test for the models. 

The models presented here are not optimal. Improved 
models could lead to a better understanding of the physics 
involved in the obscuring torus around an active nucleus. 
They could be made self-consistent with a scale height 
equal to the local radiative pressure scale height. The in- 
put spectrum could be adapted to better reflect the emis- 
sion spectrum from an accretion disk, which often shows 
a strong bump in the soft X-rays. In addition, there is ev- 
idence that star-formation can occur in the torus, which 
will increase the scale height of the disk. We plan to take 
this into account as well. Finally, we have not treated in- 
direct heating by the hot gas around the accretion disk. 
This gas absorbs the X-rays and re-emits the energy in the 
UV. This UV serves as an additional heating source for 
the dust. However, improvement of the models is only rel- 
evant when better observations are available to constrain 
the models, and when we have a better understanding of 
the contribution of star-formation to the infrared SED. 
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